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We describe the first lattice determination of the strong coupling constant with 3 flavors of dynamical quarks. 
The method follows previous analyses in using a perturbative expansion for the plaquette and Upsilon spectroscopy 
to set the scale. Using dynamical configurations from the MILC collaboration with 2+1 flavors of dynamical quarks 
we are able to avoid previous problems of haivng to extrapolate to 3 light flavors from and 2. Our result agrees 
with our previous work: af^{Mz) = 0.121(3). 



1. INTRODUCTION 

Precise determinations of the strong coupling 
constant are important for strong interaction phe- 
nomenology and searches for new physics. Lattice 
QCD provides a good methodology for such a de- 
termination [pj. 



2. IMPROVED STAGGERED QUARKS 

The staggered quark formulation allows for 
much faster simulation of dynamical quarks than 
other formulations. It maintains a remnant chi- 
ral symmetry at the cost of not entirely removing 
the doublers, which appear as additional 'tastes' 
of quarks. The presence of the doublers gives rise 
to large lattice artefacts through taste-changing 
interactions. The realisation that these can be 
systematically reduced by a Symanzik improve- 
ment programme opens the way to generating en- 
sembles with light dynamical quarks and small 
lattice artefacts. The MILC collaboration have 
generated a number of dynamical ensembles us- 
ing the Asqtad action for the dynamical quarks 



and either 2 dynamical flavors, 3 degenerate dy- 
namical flavors or 2 degenerate plus 1 heavier dy- 
namical flavors, which should be close to the real 
world In the 2-1-1 flavor case they have been 
able to reduce the light dynamical mass to 1/5 of 
the heavier one (which is approximately at nis). 

The Asqtad staggered action includes a set of 
staples in place of the simple 1-link derivative, 
to suppress gluon exchange with momentum 7r/a 
which couples the tastes together. These staples 
include a correction for low-momentum a? arte- 
facts introduced by this procedure. In addition 
there is a 3-link derivative (Naik) term to improve 
rotational discretisation errors through a?. The 
Asqtad action then has agO? rotational errors and 
a^a^ taste-symmetry breaking interactions. This 
is combined in MILC ensembles with a gluon ac- 
tion improved to the level of a^a? discretisation 
errors. 

We have calculated the T spectrum on these 
ensembles using the NRQCD action for the b 
quark 0. 
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3. NRQCD 

The physics of the bottomonium (T) system is 
that of non-relativistic b quarks [v^ jc? « 10%) 
and so the NRQCD action is appropriate 0. We 
use a Hamiltonian correct through a^w^/c* where 
all terms are tadpole-improved using an estimate 
of the Landau gauge link For the purposes 
of determining the lattice spacing to set the scale 
for (Xa we need only the radial and orbital exci- 
tation energies. These are experimentally well- 
determined and have the advantage that they are 
to a good approximation independent of heavy 
quark mass so that sophisticated tuning of the h 
mass is not necessary. In addition they are good 
probes of the dynamical quark content of the the- 
ory because the T has no valence light quarks. 

The leading contribution to the radial and or- 
bital sphttings is at ©(w^/c^) and so the expected 
systematic error from higher order (e.g. OLgV^ jc^ 
or w^/c^) terms not included in the NRQCD ac- 
tion is expected to be at the 1% level. The split- 
tings can be determined very precisely on the lat- 
tice because of the speed of NRQCD correlator 
calculations. 
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Figure 1. Q;p(8.2GeV) for various MILC ensem- 
bles with differing numbers of dynamical flavors. 
Two different quantities are used to set the scale: 
the T IP-IS splitting and the T 2S-1S splitting. 
Errors from setting the scale are included. 



4. DETERMINATION OF as 

We use the perturbative expansion of the pla- 
quette in terms of a physically well-motivated 
coupling, ay((7), defined from the heavy quark 
potential: 

y{i) = -^^vii)- (1) 

We then define a p (q) by the following expansion 
of the logarithm of the plaquette, where ap is 
equal to ay to the order (second) to which we 
are working here. 

- lii{plaq) = 3.0682(2)ap((7*)[l + (2) 
ap{-0. 770(4) - 0.09681(9)n/}] 

The numbers in brackets indicate numerical un- 
certainties in the perturbative coefficients from 
the numerical integration involved in their calcu- 
lation. This perturbative expansion is for 1-loop 
Symanzik-improved gluon action using tadpole- 
improvement with uq from the plaquette and nf 



flavors of Aqtad improved staggered quarks, as 
appropriate to the MILC ensembles. The scale, 
q* , is set using the BLM scheme and is 3.33/a 
in this case. The perturbative expansion is well- 
behaved and shows small corrections for non-zero 
n/. 

a p can be converted to aj^g using the equality 
of ap and ay and the continuum relation between 
ay and a-jj^-. 

- + 2^ + 0{al)] (3) 

5. RESULTS 

Figure 1 shows results for a p evolved to a com- 
mon (and arbitrarily chosen) physical scale of 8.2 
GeV for various MILC ensembles. Notice the 
clear separation of results for n/ = 0, 2 and 3. 

The figure shows results using either the T or- 
bital splitting or radial splitting to set the lattice 
spacing. The orbital splitting is calculated from 
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the energy difference between tfie l^Pi {hh) and 
(T) states on tfie iattice. Tlie experimentai 
result for the l^Pi mass is not known and so is 
taken as the spin-average of the masses of the l^P 
(Xfc) states. The difference is expected to be very 
small. The radial splitting is the difference in en- 
ergies of the 2^5*1 (T') and l^^i states, which is 
experimentally well-known. 

The value for ap{%.2GeV) agrees well between 
the two determinations on the 2-1-1 dynamical en- 
semble but not on the others. This reflects the 
fact that the lattice spacing will depend on the 
quantity chosen to fix it if the vacuum does not 
have dynamical content which is close enough to 
the 'real world'. 

We can take the result for ap on the 2+1 dy- 
namical ensemble, convert to ajjg and run to 

Mz- The result for a^iL is then 0.121(3). The 
dominant source of error (0.002) is the perturba- 
tive error from unknown higher orders in pertur- 
bation theory in eq. 3. This result agrees with the 
world average and with our previous results using 
dynamical (unimproved) staggered quarks 

Other sources of systematic error give small ef- 
fects. We can see from Figure 1 that the nj — 
3 results show very little dependence on dynami- 
cal light quark mass. Extrapolating ap with 
would change the final result by less than 0.001. 

The lattice spacing dependence of the result is 
expected to be very small but we currently have 
T results only on the coarser of the MILC set 
of ensembles. We can estimate the dependence 
from other results, however. Figure 2 shows the 
dependence of ap in the quenched approximation 
on for improved and unimproved glue ||^ . Here 
we have used tq to set the scale (and quote ap 
at a scale of IS/r-g) because it can be fixed very 
precisely. Its disadvantage is that it is not clear 
how to convert it to a physical result. The fig- 
ure shows that unimproved glue has very small 
dependence on a and improved glue (appropriate 
to the MILC ensembles) even less. We fold this 
lattice spacing dependence along with quark mass 
effects above and NRQCD systematics into a final 
error of 0.003. 
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Figure 2. ap in the quenched approximation at a 
scale of 15/ro for improved and unimproved glue. 



6. CONCLUSIONS 

The first calculation of as with 2-1-1 fiavours 
of dynamical quarks gives results in agreement 
with previous work, but without extrapolating in 
Uf. Further improvements, including the calcula- 
tion of the next order in perturbation theory, are 
planned to reduce the error even further. 
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